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Abstract 

The electronic properties of a Kondo impurity are investigated in a magnetic field 
using linear response theory. The distribution of electrical charge and magnetic polar- 
ization are calculated in real space. The (small) magnetic field does not change the 
charge distribution. However, it unmasks the Kondo cloud. The (equal) weight of the 
d-electron components with their magnetic moment up and down is shifted and the 
compensating s-electron clouds don't cancel any longer (a requirement for an experi- 
mental detection of the Kondo cloud). In addition to the net magnetic polarization of 
the conduction electrons an oscillating magnetic polarization with a period of half the 
Fermi wave length is observed. However, this oscillating magnetic polarization does not 
show the long range behavior of Rudermann-Kittel-Kasuya-Yosida oscillations because 
the oscillations don't extend beyond the Kondo radius. They represent an internal 
electronic structure of the Kondo impurity in a magnetic field. 

PACS: 75.20.Hr, 71.23.An, 71.27.+a 
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1 Introduction 



The properties of magnetic impurities in a metal is one of the most intensively studied prob- 
lems in solid state physics. The magnetic impurity combines two fascinating phenomena, (i) 
the formation of a magnetic moment (MM) in a host [1] , [2] due to the Coulomb interaction 
and (ii) the formation of a non-magnetic singlet state at low temperatures [3] , [1] , [5J , [6] , [7J , 
[S], [S], [ID], [EE], [12], [131 • The electronic structure of this singlet or Kondo ground state 
has been investigated using a large number of sophisticated methods, for example: scaling 
[H], renormalization [15] . [T6] . [T7] , [T8] Fermi-liquid theory [19J, [20] . slave-bosons (see for 
example [21]), large-spin limit [22], [23], local moment approach [21], [25], Bethe ansatz 
[26], [27], [28]. The author has added in recent years a new approach, the FAIR method, 
(FAIR for Friedel Artificially Inserted Resonance) which yields a very compact approximate 
ground state of the Friedel- Anderson and the Kondo impurity [21], [30], [3T] . 

One of the most controversial aspects of the Kondo ground state is the so-called Kondo 
cloud within the radius tk where tk is called the Kondo length 

hv F 

r K = — — (1) 
k B T K 

{k B T K =Ek= Kondo energy, vp = Fermi velocity of the conduction electrons). 

The idea is to divide the ground state of a Kondo impurity into two parts with 
opposite d-spins. The question is whether the conduction electrons compensate the magnetic 
moment of the d-electron forming of a Kondo cloud in each part. This question has been 
highly controversial over the last 30 years (for references see for example |32j). 

Recently the author [33] calculated the internal electronic structure of a Kondo impurity 
in real space using the FAIR solution of the Kondo ground state. This solution consists of 
only four Slater states and is well suited to calculate spatial properties which are difficult to 
obtain with other methods. The calculation showed that indeed a Kondo cloud of s-electrons 
compensates perfectly the magnetic moment of the d-electrons in the Kondo ground state. 

However, to study the polarization Kondo cloud experimentally one needs a finite mag- 
netic field because in zero magnetic field the polarization clouds of the magnetic components 
cancel. The magnetic field should be small enough so that the Kondo state is neither de- 
stroyed nor dramatically changed by the field. Therefore, in this paper I will use linear 
response theory to calculate the electronic structure of a Kondo impurity in a small mag- 
netic field. 

The magnetic properties, such as susceptibility, of the Kondo ground state in a magnetic 
field have been studied theoretically in a number of papers [31], [35], [3S], [31] • However, to 
the best of my knowledge this is the first investigation of the effect of a magnetic field on the 
electronic structure in real space. The number of real space calculations in the Kondo ground 
state is still rather limited. Recently Affleck, Borda and Saleur [3H] calculated the Friedel 
oscillations in the vicinity of a Kondo impurity by means of NRG calculations (Numerical 
Renormalization Group). The author repeated their calculation with the FAIR method and 
found good agreement with the NRG approach. 
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2 Theoretical Background 
2.1 The FAIR approach 

The basic idea of the FAIR method can be best explained for a Friedel resonance with the 
Hamiltonian 

N-l N-1 

H F = J2 Sv ^° v + Edd]d + V sd{v)[d ] c v + eld] 

It has been shown by the author [3H] , [10] that the exact n-particle ground state of a Friedel 
Hamiltonian can be expressed as by the sum of two Slater states. 



^>f = [Aal + BSj ]\a\^ (2) 



n-1 

n< 

i=l 

The state = ^a^cl is a localized state which is composed of the states cj,. The states 
a] | represent an iV-dimensional orthonormal basis representing the same Hilbert space as 

^ tA oforTYii n qo nnirinolir 4* It q full Kncic J /^i ' 



a 



the basis {<4}. The state % determines uniquely the full basis 

In this basis the conduction electron Hamiltonian Hq = ^2^ =1 £vclc u takes the form 

H = Y,^Jl E i a \ a i + E O a oO-i + Y,i V fr (0 4) a i + a l a 

As one can see the structure of H is identical to the structure of the Friedel Hamiltonian. 
The state represents an artificially inserted Friedel resonance state. Therefore I call a\ a 
"Friedel Artificially Inserted Resonance" state or FAIR-state. The use of the FAIR-states 
is at the heart of my approach to the FA- and Kondo impurity problem. Therefore I call 
this approach the FAIR method. 

In the ground state \1/ f the (n — 1) lowest states a\ are occupied (starting at i = 1). The 
states Qq and S are mixed. Further details are reviewed in the appendix. 



2.2 The Kondo impurity in a magnetic field 

The Kondo Hamiltonian in zero magnetic field has the form 

where s = a/2 is half the Pauli matrix vector, S is the impurity spin 1/2. The conduction 
electron band is expressed by discrete Wilson states cj, which have a logarithmic energy scale 
(see appendix). J u y is the matrix element of the exchange interaction between the Wilson 
states cl and cj,, via the impurity. 

In previous papers the author introduced a very compact (approximate) solution for the 
Kondo Hamiltonian which consists of four Slater states. (The path to this solution is briefly 
sketched in the appendix). 
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ij) K = Ba^d\ |O aT H ) + fldfaoj l°6T° a |> + l°»T°6l) + l°&t°<4> ( 4 ) 

Because of the spin degeneracy of the electrons two FAIR states and and the corre- 
sponding bases | a l| an d j^l j are needed. The many-electron state |0 a |0^) (|ObfO a j_)) is 
defined as 

n-l 

io fflT o 6i )= n«] T n^i4ii*o) 

The FAIR states a\ a (b Qa .) and the d states dL a form a bound two-electron state with 
negative energy. 

The four Slater states are abbreviated as 

V K = EM ad + B^ da + C^ db + CV bd (5) 

(The interpretaton of ^/ ad is that the MM- up electrons use the basis { a ir} anc ^ the state is 
occupied (cf-f empty) while the MM-down electrons use the other basis an d the state <i[ 

is occupied (&L empty). Since the ansatz for the Kondo ground state consists of four Slater 
states its secular matrix has the dimension four. One may consider the four Slater states as 
four basis states of an effective Hamiltonian, the 4x4 secular matrix. The components of 
this effective Hamiltonian are the expectation values of the original Hamiltonian ([3]) between 
the Slater states. 

The two FAIR states = X^ a o c it an d = ^2^f^o c u determine uniquely the full basis 
| a Jj an d The two FAIR states aj, and b\ are numerically optimized and determine 

the secular matrix and the ground state ^k- 

Spin-flip processes are only included in the matrix elements between the FAIR states 
a\ )a ,b\ )u and the d}_ a states. In the following I will refer to this effective Hamiltonian - 
secular matrix as the " secular Hamiltonian" . 

If one applies a magnetic field B in z-direction then the energy of electrons with their 
moment pointing upwards is reduced by Eb = ^bB while those with their moment pointing 
down is increased by Eb- Since within this paper the orientation of the magnetic moment of 
the electrons is of primary interest I use the symbol f for magnetic moment up or in short 
MM-up and j for MM-down. 

Applying a magnetic field corresponds to an external Hamiltonian 

H ex = -EbY.v Wv\ ~ n„i) + (n d1 - n dl )\ 

where n va is the number operator for the conduction electron (Wilson) states. The in- 
teraction of the magnetic field with the conduction electrons yields essentially the Pauli 
susceptibility in which we are not interested at this point. Since the magnetic field is intro- 
duced in linear response its effect on the conduction electrons can be ignored. However, for 
finite magnetic fields this is not correct as will be discussed in the appendix together with a 
sketch of the four Slater states in a magnetic field. 
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Since the effect of the magnetic field on the conduction electrons can be discarded the 
effective perturbation Hamiltonian of the magnetic field has the form 

H ex = -E B (ji d ^ -n d[ ) = E B (n d[ - n dT ) 

3 Linear Response to a Magnetic Field 

The FAIR ground state (jSJ) of the Kondo Hamiltonian is an exact eigenstate of the (approx- 
imate) secular Hamiltonian (which is given in the appendix for J = 0.1). The ground state 
has the coefficients (B, B, C, C) (with B = 0.70313, C = 0.0626 for J = 0.1). 

If we apply a magnetic field B which corresponds to a magnetic energy E B then we obtain 
a polarization in our system. The perturbation Hamiltonian in the Heisenberg picture is 

H% (0 = E B e iHt '' h (n dl - n d1 ) e~ iHt '' h 

If one switches on the field at the time to then this alters the wave function by adding a 
perturbation state 



\Mk) = - % -e- im ' h f H% it') dt' \* K (0)) 



where (0)) is the ground state in zero magnetic field in the Schroedinger picture. 
Then we obtain 



m K ) = --e- iHt ' h E B f dt'e iHt '/ h 
h Jto 



Ba^d\ + Cd^j |0 aT H > 
+ (cb^d\+Bd\al l ) |0 6T ai ) 



The term in the square brackets yields 



(5aJ T 4 - Cd]4^ |0 aT 6i ) + (Cbl r d\ - Bd\al^ |0 bT ai > 



e -iE t'/h 



This state is perpendicular to \^k) and can be expanded in terms of the other three solutions 
of the 4x4 secular matrix 

where oscillates with the frequency E\/h. This yields for 

f dt'e mt 'l h [...] e-^' h = f ■ , - a x |^ A ) c -*<*>-*a)*/a 

Jto {^o — ti\ + 

Here a factor e vt with ry = + is inserted for convergence reasons. Then the lower integration 
limit yields zero for t — > — oo. (Afterwards r\ is discarded). 
Then one obtains for the linear response wave function 

I*** (*)) = -E b Y^A-ol x |Va> e~ lEot/h 

where E x ,o = — (E — E\) > 
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3.1 Magnetic moment 

Next we calculate the magnetic moment (in units of The operator in the Schroedinger 
picture is (n^ — n^) The expectation value of the moment (//) is 



(fj) = 2 Re (6* K (t) |(n dT - n di )\ ^ K (t)> 



We calculate 



{6# K {t)\(n d i-n dl )\V K (t)> = 
ELi^^Ae-^ |(n dT - n rfi )| * K (0) e~^^ 

-&A,0 

ELiJ^aV'aI [ELi«a |^a>] 

While the four Slater states in equ. fll]) are not orthogonal with respect to each other, the 
eigenstates |^) of the secular Hamiltonian are. Therefore one obtains for the magnetic 
moment 

(^) = 2^Ea=i^I«a| 2 (6) 

For the example with the exchange interaction strength J = 0.1 I have collected the 
secular Hamiltonian, its eigenvectors and its eigenvalues in the appendix. One obtains for 
the linear response wave function \5^k (0)) = 0.984 \ipi) + 0.158 |-0 3 ). Because of the energy 
denominators (E\ — E = 2.38 x 10 -4 , E 3 — E = 0.522) the contribution of ^3) to the 
magnetic moment is almost 10~ 5 smaller and one obtains for in good approximation 

(u) = 2Eb~^— lail 2 = (°- 984 ) 2E B = 8140 * E B 
XP/ E lfi 1 1 2.38 x 10~ 4 

3.2 Kondo temperature 

In numerical normalization group (NRG) calculations one uses the susceptibility to define 
the Kondo temperature as k B Tx = 1/ (4%). The susceptibility is, in linear response, \ = 
(fj) I {E B ). If one applies this method to the example with J = 0.1 one obtains 

k B T K = — ^ = 3. 07 x 10- 5 

In the past we have used the energy difference between the singlet ground state and 
the relaxed triplet state for the definition of the Kondo energy. For J = 0.1 this yields 
k B T^ = 2.36 x 10~ 5 . These two values are surprisingly close together if one considers that 
there is some arbitrariness in the definition of the Kondo temperature and that almost every 
theoretical approach to the Kondo problem has a different definition. 
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4 Net Spin Polarization 



The ground state of the Kondo impurity does not possess any spin- or magnetic moment 
polarization because of its symmetry between MM-up and down states. Therefore one has to 
separate the ground state mathematically into its two magnetic components to answer the 
question about the existence of a Kondo cloud. The author showed that there is a Kondo 
cloud in the magnetic components of the Kondo ground state. But this polarization cloud 
can not be observed in the ground state. To measure the Kondo cloud one has to disturb 
the ground state sufficiently so that the symmetry between MM-up and down is removed 
but not too much so that the Kondo state is not destroyed. A magnetic field of the right 
strength might do the job. 

In linear response the interference between the ground state ^ k and the linear response 
state 5^ls gives the change of the electron density and electron polarization due to a small 
magnetic field. The calculation of this interference is quite similar to the author's calculation 
of the Kondo cloud in the previous paper [33J. 

When one uses the Slater state basis (^ a d, ^da, ^db, ^bd) the ground state = V'o has 
the components (B, B, C, C) and the linear response state 5^ = a.\ (2Eb/ AE m ) ip\ has the 
components 



where (B', —B', —C, C) is the eigenvector of the first excited state and for J = 0.1 the 
components have the values B = 0.70313, C = 0.0626 and B' = 0.7055, Ci = 0.0635 (see 
appendix). 

If we denote the electron density between two Slater states |\I/ Q ) and as (\I/ Q \p (r)| ^p) = 
Pa,/3 ( r ) then the total change of density due to the magnetic field is 



Since the state \^ da ) is just the state \^ a ,d) with reversed spins the two states have the 
same charge density and the opposite spin density or polarization. The same applies for 
the pairs (p db , db ; Pbd,bd) > (Pad,db, Pda,bd) and {p da ,db\ pad,bd)- Therefore one realizes that there 
is no change of the charge density. On the other hand the polarization which is defined as 
Pa,p = Pv,a,f3~ Pi;a,/3 reverses sign when the spins are reversed. Therefore one has for example 
p Qj/ g = —p{3, a - This yields a total polarization of 





P (r) = ai 



8E B 
AE 01 



(BB'p ad>ad + CC'p dbidb + (BC + CB') (p ad:db + p da4b )) 
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If one divides by (p) = a 



\2E B / AE W this yields for our example (J = 0.1) 




[0.504 p adM + 4. 1 x 10 3 p db4b + 9. 2 x 10 2 (p ad4b + p da4b 



)] 



(9) 



where BB'ja x = 0.504. 
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The calculation of the charge density of the Slater states is described in detail in ref. 
??. A summary is given in the appendix. The calculation can be performed for the one, 
two or three dimensional case. Since the Kondo impurity couples only to a single angular 
momentum of the conduction electrons (for example I = 0) it represents essentially a one- 
dimensional problem. If one defines the charge density as a one- dimensional density 5q/5r 
(for example in three dimension the charge 5q between two spheres of radii r and r + 5r) 
then the calculation and the results are essentially identical in all dimensions. (The main 
difference is that the oscillatory part of the charge shows a dimensional phase shift of D^, 
see equ. ([TO]) ). 

The four Slater states of the ground state and the linear response state are composed of 
Wilson states if> v (£). Here £ is the the distance from the impurity in units of half the Fermi 
length so that £ = rj (Xp/2). The explicit form of tp v (£) (using the logarithmic energy 
scale with A = 2, see appendix) for v < (N/2 — 1) is 



The appendix gives a more detailed discussion of the calculation of density and polarization. 

I will divide the calculation into two parts, (i) averaging over the fast-oscillating contri- 
bution to calculate the overall spin polarization and (ii) focus on the fast-oscillating part 
which will yield oscillations of the spin polarization. 

4.1 Background polarization 

Since B and B' are much larger than C and C the dominant contribution comes from the 
first term in equ. (JHJ) - The prefactor of the second term is smaller by a factor 100. The 
prefactors of pdb and pbd are smaller by a factor 10. However, their contribution is in addition 
reduced because the bases for MM-up (and MM-down) are different in the two states and 
therefore the interference terms are essentially reduced by the square of the multi-scalar 



product ( II a I^o| EI "i&o )• This multi-scalar product is about 0.14, which reduces the 



contribution of pdb and pbd by an additional factor of about 1/50. Therefore these terms 
will be neglected in the calculation. For the discussed example (J = 0.1) the polarization is 
given to good accuracy by 



This state has a fast-oscillating component cos (7r£ (l — gjpr)) (i n space) and a slowly varying 
component sin (7r£ -^p+i) • The first component yields the Friedel and RKKY oscillations while 
the second part yields the background electron density which is given by 



p°A0 = IVv (Ol 2 = 2 





V (0 ~ 2 (fj) Pad,ad 



fa) (Pda ,da Pad,ad) 



s 



In Fig.l the integrated polarization is plotted in units 2 (/i) as a function of the log 2 £ 
where £ = 2r/\p is the distance from the impurity in units of half the Fermi wave length. 




Fig.l: The integrated polarization of the ground state in a magnetic 
field in units of 2 (//) where /i is the magnetization in units fiB- 

4.2 Polarization oscillations 

In the vicinity of a Kondo impurity one obtains Friedel oscillations which have the magnitude 
(in reduced units) 

PFr (0 -po= -^A p (J-j cos (2tt£ - D|) (10) 

where D is the dimension of the system, the coefficients Cp have the values C± = 1/ (27r), 
C*2 = 1/ (27r 2 ) and C3 = 1/ (47r 2 ) in one, two and three dimensions [38]. = Hvp/kBTx 
is the Kondo length. The period in £ is equal to 1. The function A p is a universal 

function which approaches the values for << 1 and 2 for 1. (I skipped the 

phase shift 5p due to potential scattering). I will use the magnitude of the Friedel oscillations 
as a scale for the polarization oscillations in a magnetic field. Therefore I define an amplitude 
A p (£/£k) so that the calculated polarization in one dimension is given by 

p(f) =-2^(£) co8 (^-I-^) < U > 

The calculated polarization oscillations are proportional to — cos (27r£). They are shown 
in Fig.4a-e for different distances from the impurity. Each time two periods are plotted. The 
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average distances in Fig.2a-e vary between £ = 7 and £ = 32767. The minima lie at integer 
values of £ and the maxima at half integer (£ = 2r/X F ). 
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Fig.2: Oscillations of 27rj9(^) (p(^) is the magnetic polarization) as a function 
of £ = 2r/\ F . Two two periods of the oscillations are shown for different 
distances from the impurity. The average distances are from (a) to (e): 
7, 63, 511, 4047 and 32767. 
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In Fig. 3 the envelope of the oscillations of 2n(, i.e. the maxima and minima are 

plotted as triangles. The full points give the average of 27r£ *p (£) as a function / = log 2 (£). 
The amplitude is given by half the difference between the maximum and minimum curve. 
The amplitude does not show the expected long-range behavior of — 1/ (27r£) (which 
should be a constant in Fig. 2). For small distances the average is close to zero and for 
distances of the order of the Kondo length the amplitude is roughly equal to the average. 
For distances larger than the Kondo length the oscillations fade away. As discussed below 
they don't extend far beyond the Kondo length. 




Fig.3: The up and down triangles show the envelope of the oscillations of 
27r£ * p (£) where p (£) is the magnetic polarization of the conduction electrons. 
In the normalized coordinate £ (distance from the impurity) the oscillation 
period is 1. Since the abscissa is given in I = log 2 £ there are 2 15 oscillations 
in the range 2 < I < 15. The amplitude of 27r£ *£>(£) does not approach a 
constant at large distances. 



5 Discussion and Conclusion 

A magnetic field alters the ground state of the Kondo impurity and yields in linear response 
* K = B(l-e) tt , d + B (1 + e) V da + C (1 + e') * db + C (1 - e') * bd 

with 
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Therefore one has e ~ e' ~ (/i) and the components with the d magnetic moment parallel to 
the magnetic field are relatively increased in their amplitude by the value of the magnetiza- 
tion (//). Therefore the part of the Kondo cloud which is due to ty da is now only partically 
compensated by ^ a ,d- The same applies to ty d b an d ^bd- Therefore it is not surprizing that 
the net polarization has the same shape as the hidden Kondo cloud which the author calcu- 
lated recently. One even expects the scaling factor of about 2 (//}. The polarization cloud 
extends to a distance of 2 14,9 which corresponds to 1.5 x lO^Xp being of the order of the 
Kondo length. Using the Kondo energy k^T^ « 3.1 x 10~ 5 one obtains for the Kondo length 
r K = hvplksTx = A_f/ (27ikBT K ) « 5.1 x 10 3 Ai? where Wilson's linear dispersion relation 
between energy and wave vector is used. 

On the other hand the position dependence of the polarization oscillation is at first rather 
surprizing. The amplitude of (Cd/£ d ) cos (7r£) has a maximum at roughly the position 
£ » 2 10,5 = 1.45 x 10 3 which is smaller than the Kondo length but corresponds roughly to 
the distance at which the Friedel oscillation began to develop. Since the Friedel oscillations 
are supressed at short distances it is not surprizing that the polarization osccillations are 
suppressed within this range. But why do they disapper at larger distances. Here one can 
consider two physical scenarios: 

• The magnetic field polarizes the Kondo impurity and generates a magnetic moment at 
the impurity. This magnetic moment scatters MM-up and down electrons differently 
and as a result one obtains RKKY oscillations which should depend as on the 
distance from the impurity. 

• The Kondo impurity modifies the electronic structure of the host within the distance 
of the Kondo impurity. The polarization oscillations are restricted to the region into 
which the Kondo cloud of the impurity extends. 

It appears that the polarization oscillations (for small magnetic fields) are not due to the 
scattering of conduction electrons by the impurity. Instead they show the structure of the 
Kondo cloud and therefore they are restricted to distances of the order of the Kondo radius. 



6 Appendix 

6.1 The FAIR Method 

The ansatz for the ground state of the Kondo Hamiltonian (equ. (H|) has been derived from 
the ground state of the Friedel- Anderson impurity. The Hamiltonian of the Friedel- Anderson 
impurity is given 

{N-1 N-1 ~\ 

^ e u cl a c va + E d d\d a + ^ V sd {y){d\c va + cl a d a ] > + Un d ^n dL 
u=0 u=0 ) 

In mean-field approximation the Hamiltonian one replaces n d+ n d ^. =>n d+ (rid-) + (n d +) n d - 
— (n d+ ) (n d J). After adjusting (n d+ ) and (n d ~) self-consistently one obtains two Friedel 
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resonance Hamiltonians with a spin-dependent energy of the d^-state: Ed, a = Ed + U (rid-a)- 
The mean-field wave function is a product of two Friedel ground states for spin up and down 

As discussed in section II the ground state of the Friedel Hamiltonian can be expressed 
as a function of two Slater states, one consisting of only conduction electrons and the other 
containing one d-electron. 

^ F = (Aal + Btf) n<4$o 

^ ' i=l 

The state is called the FAIR state. 

The product of the two Friedel ground states in the mean field approximation can be 
expanded and yields a "magnetic state" 

0a T H ) (12) 

where | a I j an d j are ^ wo (different) bases of the iV-dimensional Hilbert space and aj, &J 
are two different FAIR states. 

n— 1 n— 1 

i=i i=i 

The state (j!2p opens a wide playing field for optimization: (i) The FAIR states and fej can 
be individually optimized to minimize the energy expectation value (ground-state energy). 

Each FAIR state defines uniquely a whole basis |aj j , The coefficients A, B,C, D can 

be optimized. This yields a much better treatment of the correlation effects. The resulting 
state is denoted as the (potentially) magnetic state $ms- 

The magnetic state ^ms is used as the building block for the singlet state. ^ m s together 
with its counterpart where all spins are reversed yield two states ^ms (Tl) an d ^ms (IT) • 
The singlet ground state is then given by 

^SS = ^Ms(n)±^Ms(^) 

A new optimization of the FAIR states and the coefficients changes and b\ drastically 
increases the weight at low energies. 

The Kondo Hamiltonian represents a limit of the Friedel- Anderson Hamiltonian by sup- 
pressing zero and double occupancy of the d-states. This leads to the approximate ground 
state H] in chapter II. 

6.2 Wilson's states 

Wilson considered an s-band with a constant density of states and the Fermi energy in the 
center of the band. By measuring the energy from the Fermi level and dividing all energies 
by the Fermi energy Wilson obtained a band ranging from —1 to +1. To treat the electrons 



^ ms — 



Aa\^ Qi + Ba^d\ + Cd\b ] oi + Dd\d\ 
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close to the Fermi level at ( = as accurately as possible he divided the energy interval 
( — 1 : 0) at energies of — 1/A, — 1/A 2 , — 1/A 3 , .. i.e. ( u = — 1/A". This yields energy cells 
with the range {—1/A U : -1/A U+1 } and the width A v = ( u+1 - ( u = 1/A" +1 . During this 
paper generally the value A = 2 is chosen. (In this paper I count the energy cells and the 
Wilson states from v — 1 to N.) 

Wilson rearranged the quasi-continuous original electron states ifik (r) in such a way that 
only one state within each cell <L V had a finite interaction with the impurity. Assuming that 
the interaction of the original electron states (r) with the impurity is /^-independent this 
interacting state in €„ had the form 



(r) = Ec^fc (r)/V z » 

where Z v is the total number of states (r) in the cell [Z v = Z (Cu+i — (u) /2, Z is the 
total number of states in the band). There are [Z v — 1) additional linear combinations of the 
states ifk in the cell <L V but they have zero interaction with the impurity and were ignored 
by Wilson as they are within this paper. 

The interaction strength of the original basis states ifik (r) with the d-impurity is assumed 
to be a constant. Then the exchange interaction between two Wilson states ip v (r) and ip^ (r) 
is given = ^ZjZ^fZjZJ = v^C+i - /2 ^(C+i - Q /2J. 

6.3 The Slater states in a finite magnetic field 

In Fig. 4 the energy bands of the conduction electrons a ia and bj a are drawn together with 
the FAIR states al a and 6q ct and the d states S c . We have four Slater states. The first 
Slater state is occupied with (n — 1) electrons of the basis and (n — 1) electrons of 

the basis |^| jthe FAIR state and the d-state d^. In the absence of a magnetic field 
the second Slater state is obtained by reversing all moments. So in the second Slater state 
the MM- up electrons would be filled into the j basis and the MM-down electrons into 

the | a !j.| basis. However, since the magnetic field shifts MM-up and MM-down electrons in 

different directions, the basis are slightly altered to |^|| and { a i|}- 

The energies of a\ and b\ are non-zero and opposite equal. In the first two Slater states 
the al state with negative energy is occupied while in the last two Slater states the bl state 
with positive energy is occupied. The occupation probability of the last two is smaller by 
roughly a factor 100 than for the first two. 
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{a} {b} Bases {b'} {a'} 




d 



■o 




{b'} {a'} 



{a} {b} 



Fig. 4: A drawing of the four Slater states in the Kondo ground state. In each 
Slater state a FAIR state (aj or 6q) and a d-state with opposite spins 
(magnetic moment) are occupied. The energies of the FAIR states are 
non-zero and opposite equal. The brackets below the bands show which 
basis the Wilson states occupy. A magnetic field has shifted the moment 
up and down bands in different directions. 



The magnetic field moves the MM-up and the MM-down bands in opposite direction by 
±Eb- The Fermi levels of MM-up and MM-down electrons readjust by electron spin flips to 
the same height. As a consequence the Fermi momenta are slightly different for MM-up and 
the MM-down electrons. 
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In the present calculation the slight shift of the bands is neglected because the shift is 
very small compared to the band energy. The Wilson states in the shifted bands measure 
their energies from the new Fermi level. Then the energy of the Wilson states is not changed 
except for the first and last Wilson state. The small change for v = 1 which extends from 
— 1 ± Eg to — | and for v = N which extends from +| to 1 ± Eb is neglected. Since the 
FAIR bases are constructed from the Wilson states their energies are, in first approximation 
not effected by the magnetic field. 

A finite magnetic field has, however, an important effect on the matrix elements. If one 
considers two Wilson states cL and cj,, with the same energy but different MM then their 
orbital wave functions are not identical. As a consequence the scalar product between the 
orbital parts of cL and cj,, is no longer equal to one but reduced to 1 — 2Eb/ (Ch-i — CO 
(for 2Eb < (Cv+i ~ CO)- I n the FAIR method such matrix elements occur only between the 
FAIR states. Since the present calculation is performed in linear response this problem does 
not arise. But for finite fields it will be important. 



6.4 The Secular Hamiltonian and its Eigenstates 



The FAIR solution for the Kondo ground state is given by the ansatz (J5J) which consists 
of four Slater states. With this state one can calculate the expectation value of the 
Hamiltonian Hk (equ. ([3])). The minimization of this expectation value with respect to the 
four coefficients (which at this stage are considered as independent) yields a 4 x 4 secular 
matrix. In parallel the variation and optimization of the FAIR states and 6q yields the 

optimal bases an d j^l} with the lowest ground-state energy. (The calculation and 

optimization of the energy expectation value has been described in ref. |31j). 

The secular matrix possesses four eigenvectors with four energy eigenvalues. The co- 
efficients of the eigenvector with the lowest energy are the coefficients of the ground state 
IV'o) — \^k)- The other eigenvectors represent excited states l^i),]^) , IV^) which are needed 
in the linear response calculation. With respect to these eigenvectors and their energy the 
secular matrix plays the role of a Hamiltonian, which I baptized the secular Hamiltonian. 
For J = 0.1 and a Wilson basis with A = 2 and N = 50 this secular Hamiltonian is given 
below by equ. f|T3|) . 



( -3.0013 -.00019237 -.046265 -.057213 \ 

-.00019237 -3.0013 -.057213 -.046265 

-.046265 -.057213 -2.4887 -.0029434 

\ -.057213 -.046265 -.0029434 -2.4887 / 



(13) 



The four eigenvectors of (1131) are given by the four columns of the 4 x 4-matrix (1141) 
below. 

/ .70313 .70547 -.076084 .049993 \ 
.70313 -.70547 -.076084 -.049993 
.06256 .063539 .70447 -.70438 
\ .06256 -.063539 .70447 .70438 / 



(14) 
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The eigenvalues are (E , E u E 2 , E 3 ) = (-3.0056,-3.0053,-2.4856,-2.4831). The energy 
difference (E 1 — E ) is E lf0 = 2.38 x 10~ 4 . The band energies for the occupied states in basis 

|a}| is E bd (a) = -1.4988 and for the basis {&}} is E bd (b) = -1.22305. 

6.5 The wave function of Wilson's states in real space 

For the discussion of the wave functions in real space we assume a linear dispersion relation 
between energy and momentum, a constant density of states and a constant amplitude at 
r = 0. Then the interaction between the original basis states tpk (r) and the d-impurity will 
be constant for all k-states. These assumptions are the same as in Wilson's treatment of the 
Kondo impurity. We define the wave functions ipk (r) in such a way that the results apply 
for the impurity problem in one, two and three dimensions. 

6.5.1 One-dimensional case 

Let us start with the one-dimensional problem. Here we have the impurity at the position 
zero and the conduction electrons are located in the range between and L. The wave 
functions ipk ( r ) have the form ip^ (r) = y/2/L cos (kr). There is another set of eigenstates 
V^fc ( r ) — A/2/Lsin (kr). These states don't interact with the impurity at the origin. There- 
fore they don't have any bearing on the impurity problem. 

6.5.2 Three-dimensional case 

In three dimensions the free electron states can be expressed as ip^ (r) oc Y/ 71 (9, <f)) j t (kr) 
where Y™ is a spherical harmonics and I, m are the angular momentum and magnetic quan- 
tum numbers, ji (kr) is a spherical Bessel function. Its long range behavior is given by 
(1/kr) sin (kr — In/ 2). Only the states with the same / as the impurity couple to the impu- 
rity. All the other states for different / belong to the group of inert states (r). 

If one calculates the density of the wave function, integrating in the three-dimensional 
case over the spherical surface 4irr 2 and averaging over short range (Friedel) oscillations then 
one obtains in the 



• one-dimensional case: (2/L) cos 2 (kr) — 1/L 

• three-dimensional case: pk (r) = (2/L) sin 2 (kr — In/ 2) = 1/L 

In both cases one obtains essentially the same density. Therefore it is sufficient to use 
the one- dimensional approach for calculating the density of a Kondo cloud. It is equivalent 
to the 3-dimensional case integrated over the spherical surface. 
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6.5.3 The wave functions in one dimension 



While the energy is measured in units of the Fermi energy the momentum will be measured 
in units of the Fermi wave number. We assume a linear dispersion relation for < k < 2 
with 

C = («-i) 

Here k — k/kp is dimensionless. It is useful to measure distances also in dimensionless units. 
We define £ = \kpr = jjj^- Then £ gives the distance from the impurity in units of A^/2. 
The (almost) continuous states (p K are given as 

where L is the length of the one-dimensional box. The boundary condition cos (ttkL) = 
yields k — (A + 1/2) /L (where A is an integer. The maximum value of A is the integer of 
(2L) , since k is dimensionless then L is also dimensionless). Therefore we have Z = 2L 
states in the full band of width 2. 

To obtain the Wilson state we have to sum the states ip K (£) over all states within an 
energy cell. If the cell ranges from (£„ : C^+i) corresponding to a /t-range (1 + £ u ) < k < 
(1 + Cv+i) then we represent all the states in this energy interval by 

From Z u = L (Cu+i — (u) states we have (according to Wilson) constructed one state ip v (£) 
which couples to the impurity. Similarly one can construct (Z v — 1) additional linear com- 
binations of ip K (£) which are orthonormal and do not couple to the impurity at the origin. 

We denote these states as <p u j (£). They are as inert to the impurity as the states f K (£) and 
will be included in the quasi-vacuum. 

After integration the wave function of the state cj, has the form for v < N/2 

_ ^ ! ing^g) + 
For the logarithmic energy scale (A = 2) this yields for v < N/2 — 1 using ( u = — 1/2" 



W (0 = 2v^ Sm ^ + ^ cos ^[1 - ^ ) ) (15) 



(^2^) 



Similarly one obtains for in the positive energy range 



K) = 2V^^i^l cos (, { (l + JL) ) 

The two wave functions ipN/2-1 and tpN/2 are special because their K-range is the same 
as their neighbors ip N/2-2 and tp 'N/2+1 All four states close to the Fermi level have the same 
K-range of 2~ JV / 2 ~ 1 . One has to pay special attention to this complication. 
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6.6 Density of the Wilson states in real space 

The density of a single state ip u (£) is given by the square of the function ip u (£) in equ. (fT5 

2_ 8 8111 [ 2 ) /^(2 + Cu + Ul) 



1^(01 =77 T^r —7-2 cos: 

This density has a fast oscillating contribution which yields the Friedel oscillations. We 
first average over the fast oscillation (which has a period of the order of 1 in units of Xp/2). 
Then we obtain for the Wilson states for v < N/2 where (Cm-i _ Cu) = <prr 



sm 



(^2^) 



TO 

In the numerical calculation we will use (most of the time) N = 50 Wilson states. The 
different wave functions ip u (£) have very different spatial ranges and therefore very different 
densities, the lowest being of the order of 2~ 25 < 3 x 10 -8 . This means that it is not useful 
to calculate the density as a function of £ because this density varies over a range of 2 25 . 
Instead, we use the integrated density, integrated from to £. 

JO JO (7rf) 

2"+ 2 sin 2 (7™) , 

— du 

(nu) 

One realizes that a single integral yields the integrated density for (almost) all wave function 
Vv(0- The state V'o (^O lies roughly in the range £ < 2 2 , i.e. the integrated density 
q v (£) = J I -00 (^Ol 2 increases in this range to 90%. Therefore the states ipv(€) and 
ipN-i-v (ttO essentially lie in the range £ < 2^ +2 (in units of \p/2). For a total of iV = 50 
Wilson states the maximum range of the wave functions is roughly 2 N / 2 ~ 1 = 2 26 . 

We may define as a ruler a linear array / (s) where s is an integer, (—N/2 < s < N) as 

r2 a ■ 2 1 \ 

r » / sm (nu) , 
I(s) = 2 \ 2 > du 
Jo (nit) 

Then the integrated density of the state ip v in the range from to 2 l is given by 
<£(2')= r llMO 1 2 ^' = 2 f^^^du = I(l-v -2) 

JO JO (7TWJ 

I (I — v — 2) gives the integrated density of the wave function -?/v within the radius 2 l . 
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6.6.1 Interference terms in the density 



The Wilson states ip u (£) or cj, represent the free electron states in the impurity problem. 
With the impurity we express the ground state in terms of new states a\ = Z~2u=o a i c t- Their 
integrated density is given by 



./0 



The quadratic terms can be evaluated with the same ruler / (s) as before. But this 
time one has in addition interference terms ip v (£) Vv+a (0 ■ These terms depend on two 
parameters, v and A. So one needs for each A a different ruler. Furthermore the interference 
terms depend on the sub-bands of ip v (£) and ip u+ \ (£). If both states lie either in the negative 
energy sub-band (is, v + A < AT/2) or in the positive sub-band (z/, v — A > AT/2) then one 
obtains one set of rulers I (z/, A) and if they lie in opposite sub-bands then one obtains 
another set of rulers I\ (v>, A). As an example one obtains 



J ( S ,A) = 2v / 2^ / V / V ^ cos 3m 1-- 



For 7i (z/, A) one has to replace the minus sign in the cosine function by a plus sign. Further- 
more one has to treat the terms where u + \ = N/2 — 1 separately because one state lies at 
the Fermi level and has a different cell width. 



6.6.2 The net integrated density 

n-l 

If we occupy all Wilson states below the Fermi level then we obtain FJ cJ,<3>o with n = N/2 

and $o the vacuum state. This state is not really the free electron ground state. To obtain 
the latter we have also to occupy the states tp K (£) and ip K (£). They don't interact with the 
impurity but they are occupied. Therefore we define as quasi-vacuum <3>q the state in which 

- rt— 1 

all states (p K (£) and ip K (£) with k < 1 are occupied. Then the ground state is \l/ — Y\ c t^'o- 
This state has a constant electron density in real space. 

n-l 

In the presence of the impurity the new ground state ^ new = Yl a l®b m ust also contain 

i=0 

this quasi vacuum, i.e., the non- interacting states must be occupied up to the Fermi level. 
Since the inert states are occupied in an d $ n e«) they cancel out when one calculates 
the change in the electron density. The net density of the new state ^ new is the difference 
between p ($ new ) and p (^o)- Since the inert states cancel out one can ignore their existence 
during this calculation. 
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